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Abstract 

We extend the description of flavour that exists in the Randall-Sundrum (RS) model to the soft wall 
(SW) model in which the IR brane is removed and the Higgs is free to propagate in the bulk. It is 
demonstrated that, like the RS model, one can generate the hierarchy of fermion masses by localising 
the fermions at different locations throughout the space. However, there are two significant differences. 
Firstly the possible fermion masses scale down, from the electroweak scale, less steeply than in the RS 
model and secondly there now exists a minimum fermion mass for fermions sitting towards the UV 
brane. With a quadratic Higgs VEV, this minimum mass is about fifteen orders of magnitude lower than 
the electroweak scale. We derive the gauge propagator and despite the KK masses scaling as m;; ~ n, 
it is demonstrated that the coefficients of four fermion operators are not divergent at tree level. FCNC's 
amongst kaons and leptons are considered and compared to calculations in the RS model, with a brane 
localised Higgs and equivalent levels of tuning. It is found that since the gauge fermion couplings 
are slightly more universal and the SM fermions typically sit slightly further towards the UV brane, the 
contributions to observables such as e# and Am^, from the exchange of KK gauge fields, are significantly 
reduced. 
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1 Introduction 

The soft wall (SW) model was originally proposed as a potential dual to a field theory with linear confine- 
ment [1], the idea being that under the AdS/CFT correspondence bound states in the, now broken, conformal 
field theory are conjectured to be dual to Kaluza- Klein (KK) modes on the gravity side 0. When, as in 
QCD, confinement is modelled such that the quark energy grows linearly with the quark separation then the 
mass of the nth meson will be subject to 'Regge scaling', m\ ~ n. However, on the gravity side, when a 
discrete KK spectrum is obtained by imposing two hard cut offs on the space, for example by cutting off the 
space with two branes, then the KK masses will generically scale as m„ ~ n 2 J3]|4). 

In the SW model, the hard cut off in the IR of the theory is replaced by a smooth cut off. To be specific, 
one considers an AdSs space which is cut off in the UV by a brane but extends to infinity in the IR in an 
analogous fashion to the the RSII model [5 j. However, in the IR the space has a cut off due to a smoothly 
decaying background value of a 'dilaton' given by 



Although we refer to this field as a dilaton, it should be stressed this is very much a 'bottom up' model. 
When v > 1 this smooth cut off gives rise to a discrete (bosonic) KK spectrum with (6) 13 



Hence with v - 2 one obtains a Regge scaling in the KK spectrum. For the majority of this paper we 
shall focus on this scenario, although it is worth pointing out that when v < 1 one obtains a continuous 
KK spectrum, or in other words the propagator has no poles. Such a scenario is conjectured to be dual to 
unparticles. When v = 1 one obtains a continuous KK spectrum above a mass gap. 

In recent years there has been significant interest in the phenomenological implications of AdSs in the 
context of the Randall and Sundrum (RS) model 00. This is principally because the model appears to 
offer a geometrical explanation of the two apparent hierarchies of the standard model (SM). In particular 
by localising the Higgs on the IR brane, the electroweak scale is gravitionally redshifted down from the 
Planck scale and hence the model offers a potential resolution to the gauge hierarchy problem. Likewise, by 
allowing the SM fermions and gauge fields to propagate in the bulk, one can use a bulk mass term to localise 
the fermions towards the UV and away from the Higgs. Such a scenario not only offers an explanation of 
the hierarchy in the fermion masses (9l [TUJ [HI but it also offers a suppression of potential flavour changing 
neutral currents (FCNC's) lfl2l[T3ll . While this provides a compelling description of flavour it is found that a 
problem arises when one considers the heavier quarks. With the Higgs localised on the IR brane the possible 
fermion zero mode masses scale down from the electroweak scale and hence the heavier quarks typically 
have to sit quite far towards the IR brane and in a region where FCNC's are no longer suppressed. This 
results in constraints, in particular from K° - K° mixing, forcing the mass of the first KK gauge mode to be 
;>20-30TeVfm[T5][El|. 

From the holographic perspective the RS model is conjectured to be dual to a theory which is closely related 
to technicolor ifTTl [T8l [T9l . albeit technicolor without linear confinement. So a natural question one may 
ask is which of the phenomenological features of the RS model be translated across to the SW model. For 
example, can the space be stabilised such that the gauge hierarchy problem be resolved? When one stabilises 
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such spaces using the Goldberger-Wise mechanism KOll one finds an IR curvature singularity |2"T1I22"1|2"51 . 
Hence with no IR cut off, the SW model suffers from a naked singularity. In any valid solution the boundary 
terms in the equations of motions must vanish at this singularity Il24l . This provides an additional constraint, 
on the SW model, which forbids the space from being cut off too sharply. In particular it has been argued 
that a warp factor, sufficient to resolve the gauge hierarchy problem, can be generated when v < 1 E51 or 
y < 2 [24]. None the less it is interesting to note that Regge scaling, in the KK masses, occurs in a class of 
SW solutions at the limit of what can still potentially offer a resolution of the hierarchy problem. 

Of course generating a large hierarchy is only one aspect of resolving the gauge hierarchy problem. One 
must also demonstrate that electroweak symmetry can be broken such that corrections to electroweak ob- 
servables do not force the KK scale to a level at which one reintroduces a so called 'little hierarchy problem'. 
In the SW model, with the absence of an IR brane, the Higgs must necessarily propagate in the bulk. This 
gives rise to a z dependent Higgs VEV, h(z) ~ hoz". If such a Higgs VEV is too flat then it is found that 
one requires fine tuning in order to generate the correct electroweak scale, in particular it is found that to 
avoid this tuning a > 2 Il26l . However a bulk Higgs also leads to a suppression in the couplings between 
KK gauge fields and the Higgs, which in turn gives rise to a significant reduction in the constraints coming 
from electroweak observables (6l|26j|27l|28|. The extent to which electroweak constraints are suppressed is 
sensitive to how flat the Higgs VEV is although, with a nearly quadratic Higgs VEV, it is possible to reduce 
the constraints to Mkk <: 1 TeV without a custodial symmetry [26]. 

The focus of this paper is to examine whether or not the RS description of flavour can be translated to 
models with a SW: Can the model generate a hierarchy of fermion masses without a hierarchy in the Yukawa 
couplings, and to what extent are FCNC's suppressed? Before one can carry out a meaningful study there 
are a number of issues that must be addressed. Firstly, as explained in G91I301I3T1 . when the fermions are 
propagating in the bulk, the smooth cut off provided by the dilaton is not sufficient to give rise to a discrete 
KK spectrum. In order to ensure the latter, one must couple the fermions to something that is gaining a z 
dependent VEV. A natural candidate for this (used in E91 l30l l32lB is the Higgs. However here, following 
||3T1 . we introduce a generic z dependent mass term that could arise from the Yukawa couplings to the Higgs 
or some other new physics (such as couplings to the dilaton or a Goldberger-Wise scalar [33]). It should be 
stressed that such a z-dependent term must be present in any realistic soft wall model with bulk fermions. 
We then proceed to demonstrate that a hierarchy of fermion masses can be generated when a > 2. 

FCNC's occur in such a model due to the tree level exchange of KK gauge fields. The standard approach 
to computing the size of such FCNC's is to integrate out the KK modes in order to obtain the coefficients 
of four fermion operators. The computation of such coefficients will include a sum over the KK modes. 
Naively one may be concerned that, with n? n ~ n, such a sum would be logarithmically divergent, i.e. such 
processes would be dominated by the higher KK modes. It is generally thought, however, that models such 
as the present one are effective theories and unknown UV dynamics will cut off any such divergence. Either 
way, in section [3j the gauge propagator is computed and it is demonstrated that, provided the fermion bulk 
mass is sufficiently z-dependent, this divergence does not occur at tree level. 

In section @] we compute the tree level corrections to certain flavour violating lepton decays. This study 
allows us to check the validity of our approximations against earlier work [32 ] and examine some of the 
central physics. Finally, in section [5j we apply our results to a study of the quark sector. It should be stressed 
that this can not be considered a complete study. Firstly we do not consider SW models in generality but 
focus on a specific model with a quadratic Higgs VEV (a = 2) and a quadratic dilaton (i.e. Regge scaling in 
the KK spectrum, v = 2). This model has been chosen since it is believed to be close to the optimal resolution 
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of the gauge hierarchy problem, i.e. it is still able to generate the large hierarchy but with minimal constraints 
coming from the electroweak precision tests [26]. It also allows for a reasonable level of analytical control 
over the calculation. Secondly we only consider the tree level contribution to gauge mediated FCNC's in 
the kaon sector. As mentioned, some of the tightest constraints on the RS model come from the contribution 
to €k coming from the exchange of KK gluons. Hence here we primarily focus on the kaon sector. It is 
found that, due to the Higgs propagating in the bulk, all the SM fermion zero modes can sit further towards 
the UV brane than in the RS model and hence one finds that the contribution to €k is much smaller in the 
SW model. So although we focus on the kaon sector we anticipate that the same effect would cause an 
increased suppression in all FCNC's compared to that of the RS model, with the Higgs localised on the IR 
brane. However, as found in 041 . one would anticipate that an analogous suppression of FCNC's would be 
present in the RS model with a bulk Higgs. Although we also find a slight increase in the universality, of 
the couplings between KK gauge fields and UV localised fermions, which is caused by the removal of the 
IR brane. We conclude in section [6l 



2 Obtaining the Fermion Mass Hierarchy in the Softwall Model 

In extending the RS description of flavour to the SW model one faces a number of additional complications 
related to the backreaction of the Higgs on the fermion profiles. In the RS model, studies of flavour typically 
involve large scans over many possible Yukawa couplings lH4l[T5l[T6ll33ll3"6ll . In the SW model, the fermion 
profiles receive contributions from the Yukawa couplings and the bulk mass terms. Hence the same scans 
over parameter space would be computationally challenging. As explained in (7J |29j I3T1 . an interesting 
feature of the SW model is that one must couple the fermions to the Higgs VEV, or some z dependent mass 
quantity, in order to obtain a discrete KK spectrum and so one cannot simply neglect such backreactions. In 
this paper, following IT3T1I371 . we shall employ a bulk mass term consisting of a constant and a z-dependent 
part. Fortunately we find that, when the fermions are localised towards the UV, deformations in the fermion 
profile brought about by changes to the constant part of the (effective) bulk mass terms are far larger than 
those caused by changes in the z-dependent part. Hence in this paper we shall assume the z dependence as 
flavour-independent, giving a universal correction to the fermion profiles. The validity of this approximation 
is checked in this paper. In this section we shall clarify our assumptions, compute the fermion profiles and 
demonstrate when one can generate a fermion mass hierarchy. 



2.1 The Fermion Profile 

Let us begin by considering the SM fermions propagating in AdSs, 



where rf v = diag(-i ). The space runs from a UV brane at z - R to infinity but is dynamically cut off 

by the background value of the dilaton, 



As mentioned in the introduction here we shall focus on a quadratic dilaton that yields Regge scaling in the 
KK masses. In this model the KK scale will be given by Mkk = jr, while the 4D effective Planck mass will 
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be of the order Mpi ~ ^. Hence, in order to offer a potential resolution to the gauge hierarchy problem one 
requires a warp factor of 

Q^-~10 15 . 
R 

Here we shall not consider the stability of the space but simply impose this warp factor by hand and refer 
readers to [24] for specific realisations of this model. In order to yield a low energy chiral SM it is necessary 
to include both doublets Q¥) and singlets (T) under SU(2). If one couples the fermions, with a Yukawa 
coupling Y, to a bulk Higgs (assumed to be a S U (2) doublet) that gains a VEV h(z) then a single generation 
is described by 



fd 5 x Vge" 



+-(/Yr M v M r - /v M rr M T) - m t yy + f<wy + h.c.) 



(4) 



where T M = E^y A , is the Fiinfbein, y A are the Dirac matrices, V^/ = Dm+lom is the covariant derivative 
and um is the spin connection. The upper latin indices run over all five dimensions, while //, v run from 
to 3. Since in five dimensions fermions are not chiral, one can include a bulk mass term M. 



We decompose the fermion field into two Weyl fermions such that *F 
then make the KK decompositions 



ip L + i/r R and y 5 if/ L ,R = +^l,r and 
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requiring 



dzif (n) f (m) + f (m) \ = 6 n 

J R 



If one defines the four dimensional effective mass by the Dirac equation iy^d^rL 
equations of motion are given by Q [30j [32l 
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(7) 



It is now possible to see the double role of the Higgs. Firstly it is giving mass to the fermion zero modes 
which are associated with the SM particles. Secondly, in order to arrive at a discrete fermion KK spectrum, 
it is necessary for Q to be convergent and hence the fermion profiles must go to zero as z — > °°. Here this 
is achieved by coupling the fermions to the Higgs such that lim^oo ^ = oo [30]. One can also see the 
difficulties that would arise in scanning over many anarchic Yukawa couplings, since for each 3x3 Yukawa 
matrix one would have to solve the six coupled differential equations and normalise all the solutions. 



Hence here we shall follow [31 ] and introduce a z dependent, flavour-diagonal, bulk mass term 

s = Jd 5 x v^-* (*(;r M v M - M(zw) , 



(8) 



where 



\ c o ci z u 
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There are two possible interpretations of such a mass term. One possibility is to view such a mass term as 
arising from diagonalising ([7]), and hence one would anticipate that there would be nontrivial relationships 
between cq and My/r and between ciz a and Yh(z). This equivalence has been discussed in [30]. Alterna- 
tively one could attribute the term to some new high scale physics, in which case the exponent of the mass 
term a need not be the same as that of the Higgs VEV. Here we shall try to remain as open minded as 
possible as to the underlying physics, keeping in mind that, in order to arrive at a discrete fermion spectrum, 
c\ cannot be zero and hence any realistic scenario must give rise to such a term. Through a slight abuse of 
notation, we shall refer to the cq term as a bulk mass parameter and the c\ term as the term arising from 
the Yukawa coupling. The advantage with introducing such a mass term is one can now find analytical 
expressions for the fermion profiles. In particular after making the same KK decomposition (f5]) such that 



f 

Jk 



and the equations of motion are now 



where the ± act on and respectively. The zero mode profile (m„ = 0) can now be solved for 
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(12) 



where T{a, x) = f°° e~ 1 f~ x dt is the incomplete gamma function. Note the zero mode only exists if ^ > 
for fa (and tf/R has UV Dirichlet boundary conditions) or — < for (and fa has UV Dirichlet boundary 
conditions). Also, as in the RS case, the ^/l(iPr) profile will sit towards the UV when cq > i (co < -i). To 
find the profiles of the KK fermions one can combine (fTTT) to obtain 
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Imposing boundary conditions, obtained from requiring (fTTI) to allow a zero mode, this can be analytically 
solved when a - 2 to give ||3T1 
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(13) 



where U(a,fl,x) are confluent hypergeometric functions or Kummer functions. The KK masses of these 
fields have been included in table Q] and are in agreement with IT3TH . 



2.2 The Fermion Masses 



Having obtained the fermion zero mode profile we will treat the fermion mass, arising from the coupling to 
the Higgs, as a perturbation. I.e. we let the fermion zero mode masses be approximated by 

r°° r 

Mij= dz-Yijh&flWiiz), (14) 

JR Z 
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where i, j are flavour indices and f LR is the zero mode profile (Hzl ) with cq = c Q ' and ci = Cj' . Here we 

parameterise the Yukawa couplings by Fy = Ajj V/?, where /ly are taken to be complex and order one. The 
Higgs VEV is assumed to be of the form _7] 

h(z) = hoR-t (jpf , (15) 

where ho is a dimensionless constant. If one assumes that the ci term in __) has arisen from the Higgs VEV 
then there will be a non trivial relationship between the c { ' values and hoA/j. Since the only way to know 
such a relationship would be to solve __>, one would wonder whether the situation has been improved. How- 
ever fortunately it is found that, when the fermions are sitting towards the UV brane, the four dimensional 
effective couplings, that govern the low energy phenomenology, are far more sensitive to changes in the 
Co bulk mass parameters than changes in the ci parameters. This allows us to make the approximation of 
assuming a universal c\ value, i.e. assuming that c| ! = -c^' = c\. 

Turning back to the fermion masses (fT4l and substituting in (fT2l gives 



Mu = h ^ ! -f-- V (16) 



l=a, Li C ° l\-2c Li 2c Li \ ( l+2c Ri -2c R ' 

A( \ a ' a I \ a ' a 



Li , Rj 



Note that the integral in (fT4l) is convergent only when a > 0, (c^' - c^) > and c ° a rL ° > 0. One can 
also see that the Co terms always appear in the exponent while the c\ terms appear in the base. If one now 
assumes universal c\ values then this expression simplifies to 
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'R'\2 Cl ) f— L 
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+L oJ 2c 



M iJ = h 0-£\ — ) j (17) 

_fo_ 2MQ- a p / _____ f_LQ- a 
a ' a / I or ' a / 

This should be compared with the analogous expression for the RS model with the Higgs localised on the 
IR brane and /. ( , 0) p , = Ne TCL - Rkr ifTTlFi 



ij = ^—f L (R)fl(R) = 7(1 - 2<)(1 + 2c j R ) ; a (18) 



-y(Q 1_2c i - 1) (a 1+2c « - 1) 
where v w 174 GeV is the Higgs VEV. 



The ranges of possible masses have been plotted in figure _Q One can see two significant differences from 
the RS case. The first and most obvious, is that as one localises the fermions closer and closer to the UV one 
hits a minimum fermion mass. The location of this minimum is very sensitive to both the warp factor and 
the form of the Higgs VEV. In particular one would struggle to generate the hierarchy of fermion masses 
with a linear Higgs VEV or, for that a matter, a small warp factor. What is particularly interesting is that 
the optimum solutions to the gauge hierarchy problem, found in ESI , appear to correspond to O « 10 15 and 
a ~ 2. This would put the minimum fermion mass about fifteen orders of magnitude lower than the mass of 

4 Here we have used the metric ds 2 = e~~ lkr rf 'dx^dxy - dr 2 with r e [0, R] such that Q. = e kR 
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(a) The RS Model 



(b) a = 1 




(c) a = 2 (d) a = 4 

Figure 1: The hierarchy of fermion masses for different Higgs VEV exponents. Here c\ - 1 but it is found 
that the results are not sensitive to the value of c\, while cl,r refer to the cq constant bulk mass terms. We 
take Q - 10 15 and M KK = f = 1 TeV. 

an IR localised fermion. If one took h^M^K ~ 100 GeV then this would correspond to a minimum fermion 
mass of the order 0.1 meV which is roughly the scale of the observed neutrino masses. It is curious to note 
that this approximately corresponds to models which are found to have minimal electroweak constraints 
J26J. Hence the majority of this paper will focus on models with a = 2. One could obtain a lower mass by 
considering a split fermion scenario, in which the left-handed fermions are localised towards the opposite 
side of the space to the right-handed fermions. However one would suspect that this scenario would give rise 
to large FCNC's and hence would not be phenomenologically viable. Note that this potential relationship 
between the neutrino masses and the warp factor, in scenarios with a bulk Higgs, was first pointed out in 

EE 

The second difference, from the RS model, is the gradient of the slope. In both the RS model and the 
SW model the possible masses drop away exponentially. However, in the SW model, the exponent is no 
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longer constant. This leads us to anticipate small mass hierarchies for light fermions, such as neutrinos, and 
larger mass hierarchies between heavier fermions, such as the the top and bottom quarks. This reduced cq 
dependence also has implications in the suppression of FCNC's which we shall consider in later sections. 

It should be noted that, although here we have treated the Yukawa couplings as a perturbation, these results 
are in good agreement with the corresponding plot in [30] (with a = 2 and Q. = 10 3 ). One would also expect 
that these differences would be present in any model with a bulk Higgs (see for example II331 ). Since this 
notion of a minimum in the 4D effective Yukawa coupling is generic to all couplings between fermion zero 
modes and profiles sitting towards the IR. It is, for example, well known that the same effect is present in 
the fermion couplings to KK gauge fields (see figure [2(a)] ). 



3 The Gauge Position / Momentum Space Propagator 



Having demonstrated that, with a > 2 and Q - 10 , one can generate the hierarchy of fermion masses 
we shall move on to look at FCNC's. In the RS model, the tree level exchange of KK gauge fields gives 
rise to FCNC's which provide some of the most stringent constraints on the model. The computations of 
such processes involve a sum over the KK tower (~ zZn l/ m n) which, as already discussed, are potentially 
logarithmically divergent with m% ~ n. Even if such processes are not divergent, due to the slower conver- 
gence with n, it is clearly necessary to include more KK modes than in the RS model in order to arrive at 
an accurate result, although it is not clear how many need to be included. With this in mind we shall now 
compute the 5D gauge propagator that inherently includes the full KK tower. For the remainder of this paper 
we shall restrict ourselves to studying models with a quadratic Higgs VEV (a = 2). 



3.1 The Photon Propagator 



We start by considering a U(l) gauge field in the space Q 

S=j d 5 xjge-°(~F MN F MN \, 



(19) 



where M, N run over the five space time indices and ji, v run over the four large dimensions. Integrating by 
parts and introducing a gauge fixing term of the form 

-Cgf = e~ m (df.Af - ^°|5 5 (f e-®A 5 )f, one gets 



e~® R 

— A 5 d 2 A 5 + — /is <>=, 

2 z 2 z ' 



— /\s <h \e —0% I c " — A; 



(20) 



In this paper we will work at tree level in the unitary gauge (£ — > oo) and so neglect contributions from 
ghosts and the unphysical A5 goldstone bosons. We refer the reader to [39] for a discussion of alternative 
gauges. After Fourier transforming with respect to the four large dimensions, such that = id^, the gauge 
propagator is given by 

(AW) = -iG p (z,z , )[rf v -?^P\, (21) 
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where 



(e-^-ds) + /> 2 )g„(z,z') = e*i^ - z')- (22) 



z 

The point of all this is that in order to obtain a 4D effective theory one just has to integrate out the Green's 
} p (z, z')- With a dilaton of the form CD = ^ 



2 

function G p {z, z')- With a dilaton of the form CD = -pj, the most general solution of (1221) is 



G {z z>) = i A z 2 M(\ - a,2, 4~p + BfUQ. - a,2, p?) if z < z' 
pZ ' Z \cz 2 M(l-a,2,^) + Dz 2 £/(l-a,2,^) if z>z' 



(23) 



where A, 5, C, D are constants of integration. We have also introduced the quantity 

R ,2 p 2 



a - 



4 ' 

while M(a,fi,x) and U(a,B,x) are again Kummer function^]. Throughout the remainder of this paper we 
will repeatedly use relations taken from B01 |4T1 1421 . Following 11391 , we introduce u = min(z, z') and 
v = max(z,z') and impose the continuity condition that matches the two solutions at z = z', i.e. G p (u, u) - 
G p (v, v), obtaining 

G p (u, v) = Nu 2 v 2 UmIi - a, 2, ^ j + BU (l - a, 2, ^ 

x{cM(l-a,2,^ + Du{l-a,2,^j 
Here N is a normalisation constant found by integrating over d22b to get the 'jump' condition 

/ 2 



/ z > 2 \ z i 

Kmd z G p (u,v)f/J=ex V ^—j-, 



giving 

T(l - a) 

N = 



2(BC - AD)RR' 2 ' 
where we have used the relation 



x(M(l -a,2,x)U(-a,2,x) + (1 + a)t/(l - a, 2, x)M(-a, 2, x)) 



T(l - a). 



To fix the constants of integration we impose Neumann boundary conditions on the UV brane, i.e. d u G p {u, v)\ u= r 
0, to get 

A(a) = (Q~ 2 - a)U{\ - a, 2, Q" 2 ) - U(-a, 2, QT 2 ), (24) 
B(a) = {a- QT 2 )M{\ - a, 2, Q~ 2 ) - (1 + a)M{-a, 2, Q" 2 ). (25) 

Since there is no IR brane, the determination of C and D is a little more subtle. Here we replace one of the 
boundary conditions with a 'normalisability' condition that dictates that the propagator be comprised of KK 
modes which are normalisable with respect to J R dz e~® jf 2 - 1 . This condition implies that 



5 Here we use the notation of |40|, although M(a,/3,x) can be alternatively denoted by iF[(a,fi,x) or ®(a-,yS, x), and likewise 
U(a,j3, x) can be denoted jr" 2 F (a, l+a-fi); ; - 1 ) or x). 
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Hence we require that the integral (1261) be logarithmically divergent. For large x, Kummer functions scale 

US) 

nay 



as M(a,B, x) ~ ■^le x x a @ and U(a,B, x) ~ x a . The integrand will then scale as 



dz e-®-G p (z,z) ~ J dz.UCeR"-z~ 4a ~ l + (AD + BC)z7 { + BDe r'^z 4 "' 1 I . 

The last term is clearly convergent, and so with A and B already fixed, only by setting C = can one obtain 
a logarithmic divergence. This then results in D being arbitrary and the full propagator being given by 

T(l - a)u 2 v 2 tt u 2 \ ( u 2 \\ ( v 2 \ 

GpiU > V) = ~ 2ARR' 2 [ AM [ l - a ^-^) + BlJ [ l - a > 2 ' V2)) U [ l ~ a > 2 ' rT-i) ■ W 

The KK masses will be given by the poles of the propagator, i.e. when A(a) - 0. G p (u, v) will have the same 
form for the gluons, whereas the W/Z propagator will be slightly deformed. 



3.2 The W/Z Gauge Fields 

Before considering the W/Z propagator it is useful to first consider the individual KK modes. Working post 
spontaneous symmetry breaking, we add to ( fT9l a mass term of the form 



AX = X -Kzfg 2 A^, 



where g is the 5D coupling to the Higgs (or for the Z boson jh(z) 2 (g 2 + g' 2 )A (U A /J ) and h is the Higgs VEV 
(fT3T ). Carrying out the usual KK decomposition, A^ = f n (z)A^ \xf), such that 



f 

Jr 



dz e~®—f n f m = 8 nm , (28) 



JR 

the gauge profile will then be given by 



\R' 2 ' :: /" 2A»:- U" ' '"" f " ~ °' <29) 
When a = 2 this can be solved to give 

Mz) - A^z 2 exp - 0\ ull- an, 2, ^-^J , (30) 
where we have introduced the quantities 



8 2h l , , m 2 n R' 2 

£ = \ — - + 1 and a n = — A . 

h \ 2R 4£ 

Note on sending g 2 h 2 —> 0, £ —> \ and bearing in mind that U(a,B,x) = x Y ~^U(\ + a - B, 2 - B, x), one 
regains the gauge profiles for massless gauge fields found in 0. We have also imposed the normalisability 
condition (T28T ) resulting in the coefficient in front of the M(a,B,x) part of the solution being set to zero. 
Imposing Neumann boundary conditions on the UV brane (dsf n \ l= R- 0) then gives 

((1 + £)Q7 2 - 2a n ) U(l- a n , 2, £QT 2 ) - 2U {-a n , 2, £QT 2 ) = 0. 
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Once again as £ — > 1, the KK masses are the same as those given by the poles of the propagator (T2"4"1 ). £ can 
then be found by solving 



2<T 



(7 



W '£ ,2,^fl- 2 =2t/ 



With M KK « 1 - 10 TeV, then £ z « 1.28 - 1.0028 and hence when M KK > M W /z then g 2 ~ 0(R). In other 
words, provided the KK scale is not too large, one does not have to introduce any couplings that are very 
large or very small in order to generate the correct W and Z masses. 

We can now estimate ho by comparison with electroweak observables. Ideally one should compare to all 
observables in particular the Fermi constant. However, since one of the motivations for studying this model 
was its relatively small electroweak corrections, it is reasonable to just fit the gauge couplings and ho to 
three observables. In particular if one fits to the Z mass, fine structure constant, a, and weak mixing angle, 
G w , then 

Ana = g 2 s 2 / 2 , 

where .s 2 = sin 2 8 W and c 2 = 1 - s 2 . We have also introduced the flat normalised photon gauge profile 



fo 



and Ei(x) 



r dt— . 

Jx t 



(31) 



Hence for a quadratic Higgs VEV we find 



1 2 ^ C w s w 

hn ~ 



By fitting the couplings of the five dimensional gauge field, g and g', directly to observables we are assuming 
that the ratio of the 5D couplings is the same as the ratio of the 4D effective couplings. Hence we are 
neglecting (9(M^ Z /M 2 K ) corrections to both the W and Z couplings and masses. One can check the validity 
of this approximation by numerically verifying that 

&-1 



1 



= c 2 



which gives c 2 = 0.7564, 0.7716, 0.7759 and 0.7771 for M KK = 1, 2, 4 and 10 TeV. Using these relations 
IiqMkk ~ 262 - 245 GeV for Q. = 10 15 and Mkk = 1-10 TeV. Having obtained ho one can now repeat the 
analysis of the previous section to obtain the W/Z propagator using 



to give 



Op (u,v) 



^r(l-^ 2 v 2 exp(I^(l-^) + I^(l-^))/ / (u 

\AM\ 1 - a, 2, — - 



2ARR 



r 2 
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2 fti 2 

where a = ^ and 

A(a) = (o~ 2 (l + - 2a) U (l - a, 2, £Qr 2 ) - 2U(-a, 2, £QT 2 ) , (33) 
B(a) = (2a - Q~\ 1 + O) M (l - a, 2, ££T 2 ) - 2(1 + a)M (-a, 2, £QT 2 ) . (34) 

As already mentioned, for the majority of this paper we shall focus on the quadratic Higgs VEV (or = 2). 
None the less it is worth pointing out the special case of a linear VEV. With a - 1 , d29l ) can be solved to 
give 

/ 7 2 \ g 2 h 2 m 2 R' 2 

f n =Nz 2 U\l-a n ,2,^- i \, where £ = l + and a n = — ^ — . (35) 

Once again imposing Neumann boundary conditions on the UV brane gives the condition 

2(1 - I + a n )zN ({a n - t)U (l + I - a n , 2, CT 2 ) + U - a n , 2, QT 2 )) = 0, 

mL„r' 2 

which m turn implies £ = 1 H j — and 

f Q = Nz 2 U 1^1,2,^ = ^ 

Hence in the case of a linear VEV (a = 1), the zero mode profile of the W and Z gauge fields are flat. 
This can alternatively be seen by noting that ^5/0 = satisfies (|29l ) when g 2 h 2 = 2m 2 RR' 2 . In models with 
warped extra dimensions, the deformation of the W and Z zero mode gives rise to a number of significant 
constraints, for example corrections to the Zbb vertex, rare lepton decays and corrections to electroweak 
observables, in particular the S parameter. One would anticipate that a flat W/Z profile would lead to a 
significant suppression of such constraints, although one would expect that it would also be difficult to 
generate the fermion mass hierarchy with a linear VEV. 



3.3 Convergence of Four Fermion Operators 



We are now in a position to test whether or not the coefficients of the four fermion operators, of interest to 
flavour physics, are divergent (at tree level) with respect to summing over KK number. As we shall see, in 
section [51 these coefficients are determined by the integral 



Xoo /-»oo 
dz J r &'/i/ R (z)/i/R(z)G P («,v)^ /s (z')/i fi (z'). 



(36) 



Typically this integral cannot be done analytically, although one can make a small momentum (i.e. small a) 
approximation^] 

R' 2 ( (u 2 \ 
— 1-exp — 



G„(u,v) 



2R 



\R' 2 / A(d) 



B(a) \l av ^ -> 



(37) 



6 Here we have used that M(l, 2, x) = and Taylor expanded 



1 _ 2y - Y(l - a) -ln( 



where T is a digamma function and y is the Euler constant. 



5-2T(2-o)-4y-2 ln(^ 
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n 


- 1 










n 


= 2 








c\ = 


0.5 


C\ 


= 1 


C{ = 


1.5 


Cl = 


0.5 


Cl 


- 1 


Cl = 


1.5 


CO 


m„ 
Mkk 


^eff 


m„ 
M K k 


*eff 
VI/l 


m„ 
M K k 




m n 
M K k 




m„ 
Mkk 


^eff 
V2/1 


m„ 


J'eff 

Vzi 


0.3 


1.41 


1.90 


2.00 


1.34 


2.45 


1.10 


2.00 


2.62 


2.83 


1.85 


3.46 


1.51 


0.4 


1.41 


1.88 


2.00 


1.33 


2.45 


1.09 


2.00 


2.60 


2.83 


1.84 


3.46 


1.50 


0.51 


1.43 


1.89 


2.02 


1.34 


2.47 


1.09 


2.01 


2.59 


2.84 


1.83 


3.48 


1.50 


0.6 


1.48 


1.93 


2.10 


1.37 


2.57 


1.12 


2.05 


2.62 


2.90 


1.85 


3.55 


1.51 


0.7 


1.55 


1.98 


2.19 


1.40 


2.68 


1.14 


2.10 


2.66 


2.97 


1.88 


3.63 


1.53 



Table 1 : The masses and couplings of the first two KK fermion modes. Included is the effective Yukawa 



coupling Igjj for fermion fields that have a zero mode. M KK = f = 1 TeV ' Q = 10 



_ l 



\15 



This allows (1361 ) to be approximated as 



/-»CO 

J dzz~ 2c '° exp 
Jr 



<-c\z 2 



. J 3 
/ J A C 0~2 ( 



R 



r 2 



R 



r 2 



2 A 

1 _ J 111 

2 C 0' R , 2 



J \ 



R 



i 2 



1 jZ^ 

2 C °'/?'2 



+ 



JR 



((z -» z'), (i «-> 7')) 



(38) 



The point is that this can be approximately evaluated using 



f 

Jo 



a v T{p. + v) 



2 F { l,fi + v;n + 1; 



a+/3}' 



but only when c\ + cj - 1 > and c Q , < ^ (i.e. only when the fermions are localised towards the IR) 
B3l . Alternatively one can check when the integrand blows up. Making a large x expansion of the Kummer 
functions gives 



I 



NfT( 1 - a)N ■ 



2ARR' 2 



I dz I dz 

JR Jz'=z 



l-lr 1 

z exp 



< -c\z 2 
~R 7T 



■ exp 



a-l 



R' 2 



£1-1 



/ 2-2C 1 

z exp 



-c\z 



R 



r 2 



+ 



T(l-a) r\R' 2 )\R> 2 

dz' dz((z~z), 

JR Jz=z' 



(39) 



Clearly, for large z the integrand will be dominated by the exponential and so will go to zero only when, 
once again, 

(40) 



c\ + c[ - 1 > 0. 



It is straightforward to check this empirically, for the exact integrand, and it is found to hold for every case 
checked. So, although we have not been able to prove it explicitly, we strongly suspect that the coefficients 
of four fermion operators will only be convergent when (1401) is satisfied. That is to say, the four fermion 
operator coefficients will be divergent, at leading order, if the fermions bulk mass term is not sufficiently 
z-dependent, i.e. the ci parameter is too small. 

If one interprets the ci values as arising from the Yukawa couplings then one would anticipate that c\ ~ hoAij. 
With ho ~ (250 GcV)/Mkk then naively one would expect that quite large Yukawa couplings are 
needed in order to arrive at order one c\ values. However, as was discussed in lTT4l l36l . if the Yukawa 
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couplings are too large one loses perturbative control of the theory. One can use naive dimensional analysis 
to estimate an upper bound, on the Yukawa couplings, from ensuring perturbative control over the one loop 
correction to the fermion masses. In particular one requires 



l^effP A 2 

16/r 2 



< 1, 



mt, 



where Y e g is the effective coupling between the KK fermions and the Higgs. If one requires that the theory 
is perturbative up to at least the second KK mode (A ~ 2m„) then this requires that F e ff < 2k. Before one 
can compute F e ff > one must specify the details of the Higgs sector, in particular the bulk and brane potentials. 
However if we assume that the Higgs VEV is dominated by the zero mode and so approximate the Higgs 
zero mode profile by /h(z) ~ Nh(z), such that e'^^-f^ = 1 and hence N « this would imply that 



v (n,m) 
r eff 

V2 



Ti L dz i fHfi ^ * A L dz w f ^ K 



(41) 



where the fermion profiles are given in (fT3T >. These effective couplings are shown in table Q] Hence, as in 
the RS model, if c\ ~ (9(1) then one would anticipate losing perturbative control of the theory, at a scale 
lower than the KK scale, when \Aij\ > 3. Unfortunately it is difficult to make any concrete statements, about 
what scale one loses perturbative control of the theory, since we do not know the exact relation between the 
Yukawa couplings and the c\ parameter. It is also possible that the z-dependent mass arises through some 
alternative physics other than the Higgs. Due to these uncertainties in the model here we shall not investigate 
potentially interesting effects arising from large or small c\ values. Rather here we shall consider c\ values 
for which the results are reasonably independent of the specific ci value, in particular c\ - 0.5, 1, 1.5. 



c' 
c o 


c 1 


Propagator 


5 KK modes 


10 KK modes 


20 KK modes 


0.3 


0.4 


1.7550 x 10~ 8 


1.6729 x 10~ 8 


1.7158 x 10~ 8 


1.7325 x 10~ 8 


0.3 


0.65 


-1.8710 x 10 -8 


-1.7570 x 10~ 8 


-1.8202 x 10" 8 


-1.8471 x 10~ 8 


0.6 


0.65 


6.8463 x 10~ 9 


6.1236 x 10~ 9 


6.4615 x 10~ 9 


6.6265 x 10~ 9 



Table 2: In the propagator column is the quantity \ - 77 , while in the remaining columns is the quantity 

P 1 



-77 X 5 ! 1 "- 20 -4^2- Here Q = 10 2 , c x = 1, p = 10 GeV and M KK = A = 1 TeV. 

JQ H ~ ° ~ m ~ 



In practice, for large warp factors, the integral (1361 ) is difficult to do even numerically. In order to carry out 
the scans over parameter space required to study flavour, it is easier to work with the individual KK modes. 
Eq. (l36l ) can be equated to 



j _ V" 1 8nSn 
n=0 r 



m- 



where g 



JR 



dzf LIR (z)Mz)f LIR (z), 



(42) 



and f n is the gauge profile (1301 with f = 1). The question then remains, how many KK modes should be 
summed over. In table [2] we compare this convergence for a low warp factor. Note, at momenta much lower 
than the KK scale, the propagator is dominated by the zero mode which we have subtracted off. We observe 
a reasonably good convergence. In practice we sum over the first 15 KK modes and hence would expect an 
error of the order of a few percent. 



Before moving on to look at the results it is worth briefly looking at the relative gauge fermion couplings 
plotted in figure |2(a)| The only reason why a convergence in (l36l) is possible is because the higher KK modes 
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Figure 2: The relative coupling between the fermion zero modes and the first (solid lines), second (dashed 
lines) and third (dots) gauge KK modes for the RS model (black) and the SW model with c\ = 0.5 (red), 
c\ = 1 (blue) and c\ - 1.5 (green). Here Q = 10 15 and Mkk = = 1 TeV. The lower graph focuses on the 
universality of the couplings of fermions localised towards the UV brane. 
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are increasingly weakly coupled to the fermion zero modes. One can also see that the couplings of the RS 
model do appear to fall away more rapidly than those of the SW model, although it is not really possible to 
say much when considering just 3 modes. 

The scale of constraints from flavour physics is partly determined by the difference or non-universality of 
the coupling of different flavours, particularly when they are localised towards the UV brane (co > 5). This 
has been plotted in figure [2(b)] for an arbitrary bulk mass parameter of cq = 0.65. Firstly one can see the 
source of the so called RS-GIM mechanism since, when the the fermions are sitting towards the UV brane, 
the gauge fermion coupling is approximately universal. One can also see, from figure |2(b)[ that the SW 
model has an equivalent level of universality to that of the RS model. Critically, one can also see from figure 
|2(a)} that when the fermions are localised towards the UV, differences in the couplings are dominated by 
differences in the Co bulk mass term and not the c\ term. Hence here assuming universal c\ parameters does 
not significantly change the results. However as one localises the fermions further and further towards the 
IR then one can see that the couplings becomes increasingly sensitive to the c\ parameter. 



4 The Lepton Sector 



Before considering the quark sector we shall first look at the tree level decays involving just the charged 
leptons. The advantage to this is that one need only fit to the lepton masses and this will allow us to 
demonstrate the central physics a little more clearly. In light of the current experimental status of the PMNS 
matrix one is inclined to favour configurations with large charged lepton mixings (i.e. closely spaced cq 
parameters). However here we will consider scenarios with both large and small mixings. 

We shall also consider only the zero modes of the fermions. One may be concerned that, when the full 
mass matrix MM^ is diagonalised, the mixing of the zero modes would receive contributions from the terms 
that are off diagonal with respect to KK number. In the RS model these off diagonal terms are partially 
suppressed by the orthogonality of the fermion profiles (fTOl Ifl2l . One would expect a similar effect being 
present in the SW model although it would be partially reduced due to the presence of the Higgs profile. 
As in the case of the RS model, one should be particularly cautious, about neglecting the contribution from 
the fermion KK modes, when one has matching cq values. Since this gives rise to particularly universal 
couplings which can, if the KK modes are not considered, exaggerate the suppression of FCNC's. 

The process we will consider here is the tree level decay lj — > which is given by B4l 

G 2 F m^ 

Wj kkk) = (2|Cf/ + 2\Cfj\ 2 + \Dff + \Dfj\ 2 ) , 

where 

n n n n 

and 

Ul,r are the unitary mauices that diagonalise the fermion mass matrix. This process receives contributions 
from both the exchange of KK photons and KK Z bosons but in practice we find it to be completely dom- 
inated by the deformation of the Z zero mode plotted in figure [3] We do not include the contribution from 
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1.02 




0.2 0.4 0.6 0.8 1 



Figure 3: The coupling of the Z zero mode for the RS model (black) and the SW model with c\ = 0.5 (red), 
ci = 1 (blue) and c\ = 1.5 (green). Here Q = 10 15 while M KK = 1 TeV (solid lines) and M KK = 3 TeV 
(dashed lines). 



the Higgs or KK modes of the Higgs. The current experimental bounds on these processes are B31 1461 

Br(/T -> e~e + e~) < 1.0 x 10~ 12 , 
Br(r -> fi~n + n~) < 2.1 x 10~ 8 , 
Br(r" -> e~eV) < 2.7 x 10~ 8 . 

4.1 Numerical Analysis and Results 

This slightly simplified study essentially serves three purposes. Firstly, it allows us to look at the central 
physics for a relatively simple example. Secondly, it allows us to compare the results from the present 
model, with a z-dependent mass ©, with a model in which one uses the Yukawa couplings to get a discrete 
fermion spectrum 1301 1321 . Thirdly, it allows us to test the validity of the assumption of a universal c\ value. 
This is rather crudely tested by assuming a universal ci value and then looking at the results for any c\ 
dependence. By considering just these three decays, it is reasonable to just fit to the three charged lepton 
masses. However, one still has a sizeable number of input parameters. Although we assume a universal c\ 
value, we still allow for anarchic Yukawa couplings A in (fTTT ). We shall also assume real, flavour-diagonal 
co values. 

Here we take the cl (c q L) ) as input parameters and randomly generate ten 3 x 3 complex Yukawa couplings, 
\Aij\ € [1, 3], allowing the cr values to be solved for by fitting to the lepton masses B31 . 

m e = 0.511 MeV - 106 MeV ra r - 1780 MeV. 

We then proceed to generate a further 10,000 Yukawa couplings and for each one compute the branching 
ratios and lepton masses. Inevitably most of these configurations will not give the correct masses, and so 
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we take the 100 configurations which give the most accurate masses. From these 100 configurations we 
plot the average of the branching ratios in figure |4] We then repeat this for a hundred random KK scales 
Mkk = £ [1, 10] TeV, three c\ values and five cl values; 

(A) c L = [0.710, 0.700, 0.690] (B) c L = [0.750, 0.700, 0.650] 

(C) c L = [0.601, 0.600, 0.599] (D) c L = [0.650, 0.600, 0.550] 

(E) c L = [0.560, 0.550, 0.540]. (43) 

We then proceed to compute the branching ratios for the RS model using exactly the the same method. In 
the literature there always appears to be a slight debate over what should be referred to as the KK scale. 
Here we define the KK scale to be -j^ for the SW model and ^ for the RS model. However in the interests of 
comparison, when plotting the RS points, we rescale Mkk by a factor of jf^ such that the mass of the first 
KK gauge mode will, for both models, be about two times Mkk- The results are plotted in figure |4] 

Firstly, it should be noted that the (C) configuration, with c\ « 1, is in good agreement with the equivalent 
calculation done without the z dependent mass term [32]. Secondly, as one would expect, there is an in- 
creasing ci dependence as one moves the profile towards the IR, although this dependence is still small 
compared to changes in the Co parameters, e.g. going from configuration (C) to (D). None the less we find 
that this dependence is negligible for configurations (A) and (B). Bearing in mind that these configurations 
typically have values sitting further towards the UV this suggests that the assumption of universal c\ 
values is good when > 0.6 (or < -0.6). 

These figures also highlight the implications of the reduced Co dependence in the range of fermion masses 
(i.e. the gradient in figure [TJ). One can clearly see that as one localises the left-handed fermions closer and 
closer towards the UV brane (i.e. larger ci values), the scale of the branching ratios is reduced. In order 
to maintain the correct masses, the corresponding right-handed fermions must sitter closer towards the IR. 
Sooner or later this leads to a problem via either a large value of Br(r~ — > ju~// + /0 or even a difficulty 
obtaining the correct masses. However, due to this reduced gradient the corresponding Cq are not as extreme 
in the SW model as in the RS model. Hence the RS model has problematic points in parameter space 
before the SW model. For example a Br(ju" — > e~e + e~) < 10~ 18 is difficult to achieve in the RS model but 
not in the SW model. In other words the reduced Co dependence, in the SW model, gives rise to a larger 
phenomenologically viable parameter space than in the RS model. 

In a more complete study, that included processes such as ju — > ey, it would be necessary to include both 
neutrino masses as well as mixings. This would restrict the range of possible configurations (|43T > but one 
would anticipate that the basic results still hold. 



5 The Quark Sector 

As mentioned in the introduction, some of the most stringent constraints on the RS model come from 
FCNC's in the kaon sector lTT4l[T5l[T6ll47ll34l . so we wish to see how they are affected in the SW model. 
As in lfT4l[T6ll48l . we integrate out the the non-zero gluon and photon KK modes, as well as all weak gauge 
boson modes, at tree level, arriving at the 4D effective Hamiltonian [49] 

5 3 

^eff 2 = Z CiQ ' d + Z Ci ® d - (44) 
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Figure 4: Branching ratios for rare lepton decays in the RS model (solid dots) and the SW model with 
c\ - 0.5 (square), ci = 1 (star) and c\ = 1.5 (diamond). The five configurations considered are given in 
(|43T >. i.e. (A) in red, (B) in magenta, (C) in blue, (D) in cyan and (E) in green. Note the RS configuration 
(A) points have not been plotted due to the difficulty in obtaining a reasonable fit to the masses. The black 
lines indicate the experimental bound although the branching ratios for — > e~e + e~ lie well below the 
experimental bound. Q. = 10 15 . 
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Parameter 95% allowed range (GeV 2 ) 


Parameter 95% allowed range (GeV 2 ) 


Re Ci [-9.6, 9.6] x KT 13 
ReC 2 [-1.8, 1.9] x 1(T 14 
Re C 3 [-6.0, 5.6] x 10~ 14 
Re C 4 [-3.6, 3.6] x 10~ 15 
ReC 5 [-1.0, 1.0] x 10 14 


ImCi [-4.4, 2.8] x 10" 15 
ImC 2 [-5.1, 9.3] x 10~ 17 
ImC 3 [-3.1, 1.7] x 10" 16 
ImC 4 [-1.8, 0.9] x 10~ 17 
ImC 5 [-5.2, 2.8] x 10~ 17 



Table 3: Allowed ranges for the AF = 2 Wilson coefficients |[50l . 



Here 

Qf = (d L fs L )(d L7ll s L ) Qf = {d R fs R ){d R7ll s R ) 
Of = (d R s L ){d R s L ) Qf - (d L s R )(d L s R ) 

Qf = (d a R &(4s a L ) Qf = (^4X^4) 

Qf = (d R s L )(d L s R ) 
Qf = (d a R 4)(%® 

and a and yS are colour indices. Model-independent bounds on the Wilson coefficients C,, C, have been given 
in ||50"1 . as quoted in Table [3] Let us define the integral 



3 [" S*00 /~»DO 

^^(f/JW dz dz! y^z)f x {z)Gf\u,v)fi{z')fkz) 

; ;_i *JR +J R 



l iA) 

',7=1 



{u\) m KU<ry n , (45) 



where iff,x, £,cr = L,R while Ul/r are the unitary matrices that diagonalise the mass matrices and i, j, k, I, m, n 
are flavour indices. The Wilson coefficients are then given by 

C, = dtfW +^-f!-I J 2 fl (z) +— (46) 

1 g k L d L s L d L s ^ 2c 2 V2 3 7 L dL s L d L s ^ j § L d L s L d L s ' 

r - ^ii (G) + £^1t (Z) + — ir (A) iai\ 

u l g R d R s R d R s ^ 18c 2 L R d R s R d R s ^ ig^AKA ' 

Cs = ^I (G) +^4/I_I 5 2V( z ) -—i A) (49) 

*- 5 2 LdL s R d R s 3c 2 \2 3 H 7 L dL s R d R s Q L d L s R d R s ' 

where s^, = sin 2 6V, c 2 = cos 2 #w is the weak mixing angles and g s , g, e are the 5D couplings of the 
(five-dimensional) gluon, Z and photon fields. Here we are not interested in electroweak constraints, so we 
shall equate g 2 s f^ = Ana s , g 2 f£ = Ana/s^,, and e 2 f^ = Ana. 

The Ki - Ks mass difference and the indirect CP violation parameter 6k follow from the weak Hamilto- 
nian (1441) . Because the former is more sensitive to uncertain long-distance effects than the latter and because 
of the small experimental value (ejf) eX p = (2.228 ± 0.011) x 10~ 3 [45], the stronger constraints are on the 
(CP- violating) imaginary parts of the Wilson coefficients (cf Table 0). These stronger constraints typically 
translate into the strongest constraints on the model, as is the case in the RS model (or the MSSM, for that 
matter). We expect the same to hold for the SW model. The new-physics (NP) contributions to 6k and AMk 
are calculated using |[T6ll 

e K = ^ Im <K°|^= 2 |K°> , (50) 
V2(Am if ) exp eIt 
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1 TeV 


3 TeV 


10 TeV 


30 TeV 


Si 


0.407 


0.395 


0.384 


0.374 


B 4 


0.938 


0.938 


0.938 


0.938 


B 5 


-0.336 


-0.373 


-0.404 


-0.427 



Table 4: "Bag" factors Bj(fi) for the relevant hadronic matrix elements at representative renormalisation 
scales. 

AmK = 2Re<K |-H^T 2 |K°> (51) 

where ip e = 43.51° and k £ = 0.92 lISTI . The hadronic matrix elements of the four-fermion operators are 
parameterised as 

(K \Ql d (jj)\K°) 
(K°\Ql d (jj)\K ) 

<K°|<2f(//)|£ () > 

where = 497.6 MeV, fie - 156.1 MeV, and we have indicated the dependence on the renormalisation 
scale /i. Tremendous effort has gone into higher-order corrections to the SM calculation of C\, to the formula 
(1501 , and the determination of the 5,-factors, in particular the SM one, Si. In fact, some recent studies have 
hinted that the SM contribution may be slightly too small to explain the experimental value |T52ll53l . While 
this is intriguing, in view of our tree-level analysis and other approximations made, we will ignore the SM 
contributions, comparing the new-physics contribution in this model directly to the experimental values for 
6k and AMk- We also restrict ourselves to a leading-log analysis, renormalising our Wilson coefficients at 
fi = fio - 2Mkk and evolving [54] (at leading log) the bag factors given in 11551 up to [iq. Representative 
numerical values are given in Table 0] 



m Kfl 
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m K \ 2 m Kfl 

b 4 (m) , 



m d (2M KK ) + m s (2M KK ) j 4 
m K \ 2 m K f K 



m d (2M KK ) + m s (2M KK ) 12 



b 5 (m), 



5.1 Numerical Analysis and Results 



Here we wish to investigate two points. Firstly, can the SW model reproduce the correct quark masses and 
mixing angles and secondly, to what extent are FCNC's suppressed. However, even if one assumes real, 
flavour-diagonal, bulk mass parameters and universal, order unity, c\ values then one still has 18 complex 
Yukawas and 9 real c$ parameters to fit. Such a large parameter space gives rise to an under constrained 
problem or in other words the existing constraints from flavour physics can always be satisfied with sufficient 
tuning of the free parameters. Hence the relevant question, we wish to address here, is which of the two 
models requires the least tuning in order to reproduce all existing observables. Ideally one should carry out a 
full Monte Carlo analysis, although accurately carrying out the integrals is numerically too slow for this to be 
a computationally viable option. Another possible approach would be to compute the fine tuning parameter 
(introduced in [56]) including all know observables. However this approach would offer no indication as to 
the 'typical' size of a given observable in a given model. 

The approach taken here is to find points in parameter space that give the correct masses, mixing angles and 
Jarlskog invariant and then proceed to calculate the size of the additional contributions, to €k and Amx , from 
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Configuration 


L R 


d 

C R 


(A) 


c\ = 0.5 
ci = 1 
ci = 1.5 


[-0.66 ±0.04, -0.47 ±0.12, 0.46 ±0.10] 
[-0.65 ± 0.04, -0.45 ± 0.11, 0.47 ± 0.07] 
[-0.64 ±0.03, -0.43 ±0.13, 0.45 ±0.13] 


[-0.63 ± 0.01, -0.61 ± 0.01, -0.57 ± 0.01] 
[-0.62 ±0.01, -0.60 ±0.01, -0.56 ±0.01] 
[-0.62 ±0.01, -0.59 ±0.01, -0.56 ± 0.01] 


RS 


[-0.62 ± 0.01, -0.44 ± 0.05, 4.63 ± 1.98] 


[-0.60 ± 0.01, -0.58 ± 0.01, -0.55 ± 0.01] 


(B) 


ci = 0.5 
ci = 1 
ci = 1.5 


[-0.69 ±0.03, -0.52 ±0.03, 0.38 ±0.18] 
[-0.69 ±0.07, -0.46 ±0.14, 0.46 ±0.13] 
[-0.67 ± 0.01, -0.46 ± 0.12, 0.45 ± 0.15] 


[-0.66 ± 0.01, -0.61 ± 0.01, -0.60 ± 0.01] 
[-0.65 ±0.01, -0.61 ±0.01, -0.58 ±0.01] 
[-0.65 ± 0.01, -0.60 ± 0.01, -0.58 ± 0.01] 


RS 


[-0.65 ±0.01, -0.48 ±0.04, 1.09 ±0.67] 


[-0.63 ± 0.01, -0.59 ± 0.01, -0.57 ± 0.01] 


(C) 


ci = 0.5 
Cl = 1 
ci = 1.5 


[-0.72 ±0.04, -0.56 ±0.01, -0.25 ± 0.22] 
[-0.71 ± 0.01, -0.53 ± 0.08, 0.06 ± 0.33] 
[-0.70 ± 0.03, -0.52 ± 0.06, 0.30 ± 0.32] 


[-0.69 ± 0.01, -0.65 ± 0.01, -0.62 ± 0.01] 
[-0.68 ±0.01, -0.64 ±0.01, -0.61 ±0.01] 
[-0.68 ± 0.01, -0.63 ± 0.01, -0.60 ± 0.01] 


RS 


[-0.68 ±0.01, -0.53 ±0.01, -0.06 ±0.14] 


[-0.66 ±0.01, -0.62 ±0.01, -0.60 ±0.01] 


(D) 


ci = 0.5 
Cl = 1 
ci = 1.5 


[-0.74 ±0.01, -0.60 ±0.01, -0.37 ±0.08] 
[-0.73 ±0.01, -0.58 ±0.01, -0.23 ± 0.21] 
[-0.73 ±0.01, -0.56 ±0.06, -0.15 ±0.24] 


[-0.72 ± 0.01, -0.67 ± 0.01, -0.62 ± 0.01] 
[-0.71 ±0.01, -0.67 ±0.01, -0.61 ±0.01] 
[-0.71 ±0.01, -0.66 ±0.01, -0.61 ±0.01] 


RS 


[-0.72 ±0.01, -0.57 ±0.01, -0.29 ± 0.06] 


[-0.69 ± 0.01, -0.65 ± 0.01, -0.61 ± 0.01] 


(E) 


ci = 0.5 
ci = 1 
ci = 1.5 


[-0.77 ± 0.01, -0.63 ± 0.01, -0.43 ± 0.03] 
[-0.76 ±0.01, -0.62 ±0.01, -0.36 ± 0.07] 
[-0.76 ±0.01, -0.62 ±0.01, -0.21 ± 0.24] 


[-0.75 ±0.01, -0.71 ±0.01, -0.63 ±0.01] 
[-0.74 ± 0.01, -0.71 ± 0.01, -0.62 ± 0.01] 
[-0.73 ± 0.01, -0.70 ± 0.01, -0.61 ± 0.01] 


RS 


[-0.74 ±0.01, -0.61 ±0.01, -0.35 ±0.05] 


[-0.72 ± 0.01, -0.69 ± 0.01, -0.61 ± 0.01] 



Table 5: Bulk mass parameters (c# = c^') of the right-handed fermions obtained by fitting to quark 
masses and mixing angles with the left-handed fermions having the bulk configurations; (A) cl = 
[0.72, 0.64, 0.52], (B) c L = [0.69, 0.63, 0.49], (C) c L = [0.66, 0.60, 0.42], (D) c L = [0.63, 0.57, 0.34], 
(E) cl - [0.60, 0.52, 0.25]. Quoted is the mean and standard deviation taken over 50 random points, see 
text. 



the KK gauge fields. In selecting such points one should be aware of two factors. Firstly, is the point fine 
tuned, i.e. are the output observables sensitive to small changes in the input parameters. Secondly, is the 
point a particularly rare point in parameter space. In order to address the second issue, here we endeavour 
to scan over as wide a range of the parameter space as is computationally viable. For the first point we will 
compute the fine tuning parameter at each point considered. To be a little more explicit our method is as 
follows. 



• With the exception of the overall warp factor (Q = 10 ) and the KK scale, the only input parameters 
we fix by hand are the bulk mass parameters cl and the universal ci parameter. Here we consider five 
configurations of 

(A) c L = [0.72, 0.64, 0.52] (B) c L = [0.69, 0.63, 0.49] 

(C) c L = [0.66, 0.60, 0.42] (D) c L = [0.63, 0.57, 0.34] 

(E) c L = [0.60, 0.52, 0.25] (52) 

and three ci values, ci = 0.5, 1, 1.5. The five cl configurations have been chosen such that they 
give roughly the correct mixing angles. Note that for configurations with 0.74 < c^ l {c L ) < 0.60 it 
becomes increasingly difficult to get a good fit to the quark masses without including quite large bulk 
mass parameters. 
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Figure 5: The mean values of and Am^ p for the RS model (stars) and the SW model with c\ = 1.5 
(circles), c\ = 1 (squares) and ci = 0.5 (diamonds). The cl values are given in (l52l . For the SW model 
configuration (A) is plotted in dark blue, (B) is plotted in light blue, (C) is plotted in cyan, (D) is plotted in 
light green and (E) is plotted in dark green. While for the RS model (A) is plotted in dark red, (B) is plotted 
in light red, (C) is plotted in orange, (D) in yellow and (E) in dark yellow. For both the RS model and the 
SW model the mass of the first gauge KK mode will be about two times Mkk- Note plotted here are the 
average values. Typically one can always find tuned points, in parameter space, that satisfy the experimental 
constraints for all configurations and all KK scales considered. Q. = 10 15 . 
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Figure 6: The percentage of points that are consistent with experimental values. The colours and points 
correspond to those used in figure [5j Due to the fact that only a 100 points were used, a large number of the 
configurations and KK scales considered yielded either one or zero points consistent with 6k (figure [6(a)] ). 
Such configurations would require a significant level of tuning in order to be consistent with experimental 
constraints. 
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• Next we find the 'natural' c u R and c d R values. By natural we mean the bulk mass parameters that give 
the correct masses and mixing angles assuming that there is no hierarchy in the Yukawa couplings. 
For this we generate ten sets of two Yukawa matrices and for each one solve for c u R d by fitting to the 
quark masses. To avoid accidentally using a fairly extreme Yukawa, the cr values used is the median 
of these ten values. The Yukawa matrices are generated such that \Aij\ e [1, 3]. As discussed in 
section [331 we avoid using large Yukawa couplings. These average cr values are given in table[5] 

• With the nine bulk mass parameters fixed, we proceed to find 100 points in parameter space, that give 
the correct quark masses, mixing angles and Jarlskog invariant by solving for the Yukawa couplings. 
The quark masses are run up from 2 GeV to the mass of the first KK gauge mode (2Mkk) using the 2 
GeV values 

m u = 2.5 MeV m c (3 GeV) = 0.986 GeV m t = 164 GeV 
m rf = 4.95 MeV m s = 96.2 MeV m b = 4.163 GeV (53) 

and the mixing angles are l45l l50l 

V us = 0.2254 ± 0.00065 V cb = 0.0408 ± 0.00045 

V ub = 0.00376 ± 0.0002 J = 2.91^ x 10~ 5 . (54) 



By randomly generating the initial guess of the solver one would anticipate that these points would be 
spread evenly over the parameter space. However one still needs to check the level of tuning required 
to obtain masses and mixing angles. Hence here we compute the fine tuning parameter ll56l 



ABG(Oi,Pj) 



Pj AO, 



Oi A Pj 



(55) 



where the observables, Oi, run over all the masses, mixing angle and Jarlskog invariant and the input 
parameters, pj, run over the all the Yukawas. In practice we vary the Yukawas over a range of 0. 1 +0. 1/ 
(i.e. Apj). The fine tuning parameter is then taken as the maximum value with respect to both input 
parameters and output parameters. Plotted in figure|7]is the mean value taken over the 100 points. 

Having found these 100 viable points in parameter space we proceed to compute the size of the 
contributions to 6k and Amj-. The mean value, taken over these 100 points, have been plotted in figure 
[5] We have also plotted, in figure [6j the percentage of the points that are consistent with experimental 
results. Again in the case of the RS model the KK mass has been scaled by a factor of -^r such that 
the mass of the first gauge KK mode would be about the same in the two models. 

This process is then repeated for 50 KK scales, randomly chosen such that Mkk £ [1, 10] TeV, for 
both the RS model and the SW model with c\ - 0.5, c\ - 1 and c\ — 1.5 and also at each configuration 
in i 



As mentioned before, it is hoped that this approach will offer an unbiased spread of points over the region of 
parameter space of relevance to generating SM fermion masses and mixing angles. Let us turn now to the 
first question, posed at the beginning of this section, can the SW model generate the correct fermion masses 
and mixings? As one would expect the answer is yes. For both the RS model and the SW model, not one 
of the points plotted in figure [7] show any sign of significant fine tuning (i.e Arc > 10). This is of course 
because of bulk configurations (l52l taken as input parameters. Had one chosen more UV or IR localised 
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Figure 7: The fine tuning parameter giving an indication of the sensitivity of the quark masses, mixing 
angles and Jarlskog invariant to variations in the Yukawa couplings. The colours and points correspond to 
those used in figures [5] 



configurations then clearly more tuning would have been required in order to get the correct masses and 
mixing angles. 

This then brings us to the second question, related to the the suppression of FCNC's. By grouping points in 
parameter space according to their configurations of bulk mass parameters (l52l . one is essentially comparing 
points which require equivalent levels of tuning in order to obtain the correct masses and mixing angles (see 
figure|7]). Even with this comparison it is still difficult to meaningfully quantify the extent to which FCNC's 
are suppressed. None the less one can see that, for all points considered, FCNC's are considerably more 
suppressed in the SW model than in the RS model with an equivalent level of tuning. For example, one can 



see from figure [6(a)] that, in the SW model with configuration (E), in order to have about 20% of the points 
consistent with €k one would require Mkk ^ 4 TeV (corresponding to a KK gluon mass of ~ 8 TeV). Where 
as in the RS model one would require Mkk <: 10 TeV (corresponding to a KK gluon mass of ~ 25 TeV). 
Alternatively one can look at the total number of points that satisfy the 6k constraint regardless of KK scale. 
In the SW for configuration (E) with c\ = 1.5, 1, 0.5 this is about 31%, 31%, 19% respectively. While for 
configuration (A) this is about 4.2%, 2.5%, 2.6%. This can be compared to the RS model which is about 
8% for configuration (E) and 0.22% for configuration (A). 

There are a number of factors contributing to this increased suppression of FCNC's. The extent to which 
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FCNC's are suppressed is largely determined, at tree level, by how universal the gauge-fermion couplings, 
that appear in (|45T ). are. As can be seen in figure |2(a)[ as one localises the fermions further and further 
towards the UV, the gauge fermion coupling becomes increasingly universal. While this is true for both the 
SW model and the RS model one can also see, from figure [2(b)] that in the SW model, with c\ > 1 that the 
gauge fermion couplings are slightly more universal than the RS model. In addition to this in the SW model, 
due to the presence of the bulk Higgs, the fermions will typically sit slightly further towards the UV than in 
the RS model (see table [5]). This effect should be combined with the reduced cq dependence in the range of 
possible fermion masses (the gradient in figure [T]). This results in avoiding extreme bulk mass parameters 
in configurations in which either the left-handed or right-handed fermions are localised quite far towards 
the UV. For example, see for configurations (A) and (B) in table [5] for an extreme example. When all 
these effects are combined one finds that, for all points in parameter space considered, FCNC's are more 
suppressed in the SW model than in the RS model. 

One of the underlying assumptions in this analysis is that of a universal c\ value. However, upon examining 
the results one can see that, despite most of the bulk mass parameters of relevance to kaon physics having 
< -0.6, there is still a very small c\ dependence. In particular smaller ci values tend to give slightly more 
UV localised fermions and slightly less universal gauge fermion couplings. These effects are quite small 
and hence negligible when compared to varying the bulk mass parameters, cq. So here we would argue that 
these results would not change significantly if one was to relax this assumption. None the less, since the 
exact origin of the c\ term has not been clearly defined it is not clear if this number can be quite large or 
small. Here we shall leave investigation of this to future work. 

It should also be stressed that this cannot be considered a complete study for a number of reasons. Firstly 
we only consider tree level gauge mediated FCNC's. One would also anticipate additional contributions to 
FCNC's arising from, for example, the dilaton or the Higgs Il36ll57ll58ll59l . The inclusion of such additional 
contributions would inevitably enlarge the parameter space and hence make a fair comparison with the RS 
model more difficult. However one would anticipate that such fields would be IR localised and hence such 
FCNC's should also be suppressed. Secondly we also only focus on kaon physics. One can see from table [5] 
that a study involving top physics or B physics would involve fermions sitting further towards the IR where 
the assumption of universal ci values is arguably less valid. So here we will leave a more comprehensive 
study to future work. 

It is important to realise that the central physics involved in this result is not necessarily related to the soft 
wall but rather the change in the relationship between fermion masses and their positions. In particular the 
reduction of the gradient in figure [T] One would anticipate such an effect showing up, to some extent, in 
any bulk Higgs scenario. For example, a similar reduction, in the constraints from 6k, was found in the 
RS model with a bulk Higgs in [34]. However soft wall models offer a framework in which the Higgs can 
propagate in the bulk and the gauge hierarchy problem can still potentially be resolved. 



6 Discusion and Conclusions 



The primary motivation of this paper was to demonstrate that the description of flavour that exists in the RS 
model can be transferred across to the SW model. Here we believe that we have done this but there is one 
major difference. Whereas in the RS model, with the Higgs on the brane, the fermion masses are determined 
by the end points of the fermion profiles, in the SW model they are determined by the overlap integral of 
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the Higgs VEV and the fermion profiles. This gives rise to two significant differences between the SW 
model and the RS model. Firstly, although the range of possible masses still falls exponentially, in the SW 
model the exponent is no longer constant. Secondly, with both left-handed and right-handed fermions sitting 
towards the UV one obtains a minimum possible fermion mass. 

Both the gradient in the fermion masses and the minimum fermion mass are very sensitive to both the warp 
factor and the form of the Higgs VEV. One would also suspect that it would be sensitive to the form of 
the dilaton background value as well. Although here we have fixed v = 2 in order to arrive at a Regge 
scaling in the KK masses. It should also be noted that one cannot obtain the hierarchy of fermion masses 
if the warp factor is too small or the Higgs VEV is too flat. However, with a warp factor of Q = 10 15 the 
case of a quadratic Higgs VEV is particularly interesting since it results in a minimum fermion mass being 
approximately fifteen orders of magnitude lower than the EW scale. This appears to be in rough agreement 
with the observed range of neutrino masses and hence the majority of this paper has focused on this case. 
Further still it was found in [26] that an approximately quadratic Higgs VEV gave the lowest EW constraints. 

Having restricted the study to that of a quadratic Higgs VEV we proceeded to investigate the implications 
for flavour physics. Studies of flavour in the SW model suffer from a number of additional complications 
to that of the RS model. Firstly, one must necessarily ensure that the bulk fermions gain a z-dependent 
mass term in order to obtain a discrete KK spectrum. The most natural source of such a mass term is the 
Yukawa couplings to the Higgs. However here, following Bill , we use a generic bulk mass (© term which 
could or could not be related to the Yukawa couplings. The advantage to this bulk mass term is that it 
allows for analytical expressions to be obtained for the fermion profiles and hence simplifies the study. The 
disadvantage is that one does not know the precise relationship between the c\ parameters and the Yukawa 
couplings. Fortunately it is found that when the fermions are localised towards the UV and one assumes 
c\ ~ 0(1) then the results are relatively independent of c\. Hence for the physics considered here, notably 
kaon physics and lepton decays, one can assume a universal c\ value and still arrive at a reasonably reliable 
result. 

The situation receives an additional complication when one tries to estimate at what scale one loses per- 
turbative control of the calculation. With the KK masses scaling as ~ n, one would naively expect a 
tree level sum over the KK tower to be logarithmically divergent. I.e. tree level processes are dominated 
by the higher KK modes and hence one must impose a cut off in the KK number. By computing the five 
dimensional gauge propagator, it is demonstrated that this divergence does not occur provided the fermions 
gain a sufficiently z-dependent mass term, i.e. ci > 0.5. If one assumes that the Yukawa couplings are the 
source of the c\ term and that c\ ~ Aho then this would suggest that one needs large Yukawa couplings to 
avoid such a tree level divergence. However this would result in a loss of perturbative control, at next to 
leading order, at a scale much lower than the KK scale. Alternatively one could look for an another source 
of the z-dependent mass term such as couplings to the dilaton or a Goldberger-Wise scalar ll33Tl . 

None the less, assuming c\ ~ 0(1) and a quadratic Higgs VEV, we are then able to investigate the phe- 
nomenological implications of the change in the range of possible fermion masses. Due to the difficulties 
in quantifying the extent to which FCNC's are suppressed here we choose to compare the SW model to 
the well studied RS model. In other words, we are using the RS model as a point of reference. Firstly, 
a relatively simple study of rare lepton decays is made primarily to check the assumption of universal c\ 
values as well as allowing for comparison with some earlier work Il32l . With ci « 1, good agreement is 
found with models in which the backreaction of the Yukawa couplings, on the fermion profiles, is included. 
Next we proceed to examine K Q - K° mixing for both the RS model and the SW model, comparing points 
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in parameter space that require equivalent levels of tuning in order obtain the correct masses, mixing angles 
and Jarlskog invariant. It is found that, in the SW model, the fermions typically sit slightly further towards 
the UV than in the RS model, with a brane localised Higgs. Also the gauge fermion couplings are slightly 
more universal. This results in FCNC's being more suppressed in SW model than in the RS model for nearly 
all points in parameter space considered. Hence here we would conclude that, despite its computational dif- 
ficulties and uncertainties, the SW model arguably offers a more appealing description of flavour than the 
RS model. This result, coupled with the result of 1126*11271 . suggests that there is a strong case to be made for 
considering models with a bulk Higgs. 
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